September, 1995 



Quantum duality and Bethe-Ansatz 
for the Hofstadter problem on hexagonal lattice 

C.-A. Piguet, D.F. Wang and C. Gruber 

Institut de Physique Theorique 
Ecole Polytechnique Federale de Lausanne 
PHB-Ecublens, CH-1015 Lausanne, Switzerland 

The Hofstadter problem is studied on the hexagonal lattice. We first establish 
a relation between the spectra for the hexagonal lattice and for its dual lattice, the 
triangular lattice. Following the idea of Faddeev and Kashaev, we then obtain the 
Bethe-Ansatz equations for this system. 

PACS number: 71.30.+h, 05.30.-d, 74.65+n, 75.10.Jm 



Typeset Using REVTEX 



1 



Systems in external magnetic field have been of considerable interest in recent years. One 
of the most fascinating properties of these systems is the integer and fractional quantum 
Hall effect [0-0] . The essential physics of the integer quantum Hall effect can be described 
by the simple Landau problem, in which a free electron moves in a two dimensional plane 
under constant magnetic field. For the fractional quantum Hall effect, it is well known that 
the electron-electron correlation gives rise to the energy gap of the system. 

Besides the systems of electrons moving continuously on the two dimensional plane in an 
external magnetic field, the problem of free electrons hopping on a two dimensional lattice 
under external magnetic field, i.e. the Hofstadter problem, has attracted a lot of attention 
PHTTf . Recently, using reflection positivity, Lieb has provided a proof for the long standing 
conjecture that for the square lattice, the magnetic flux which minimises the energy of the 
system at half-filling is exactly ir per plaquette |10| . On the other hand, in their recent work, 



Wiegmann and Zabrodin have shown that the magnetic translations can be constructed with 
the generators of the quantum group U q (sl(2)) for the Hofstadter problem ||. With this 
representation, they obtain the Bethe-Ansatz equations for the eigenvalue problem of the 
system on a two dimensional square lattice. Later, Faddeev and Kashaev were able to 
provide a generalized approach, both for the square and triangular lattices ||. 

In two dimensions, there are three lattices of special interest: the square, the triangular 
and the hexagonal lattices. In particular, the triangular lattice and the hexagonal lattice 
are dual of each other. In this paper, we study the Hofstadter problem on the hexagonal 
lattice. Following the idea of Wiegmann- Zabrodin, and Faddeev-Kashaev, we shall find the 
Bethe-Ansatz type solutions for this system. Furthermore, we establish a quantum duality 
relation between the energy spectra on the triangular and on hexagonal lattices. 

The general Hamiltonian for non interacting electrons moving on a lattice in a magnetic 
field is: 

// (i) 

where cf , C{ are the creation and annihilation operators for an electron ai site i, tij is the 



— * — * 

real hopping matrix, and corresponds to // A(x) ■ dx, with A the vector potential and the 
integral is performed on a straight line between % and j. It thus satisfies 8ij = In the 

following, we consider the special case where the electrons hop between nearest neighbours 
only. We further assume that the system is invariant under translation, but not necessarily 
invariant under rotation. 

The Hofstadter problem is a very interesting problem of theoretical physics, in particular 
due to its distinction between rational and irrational numbers. In the following, we only 
consider the rational case when the flux per elementary cell of the lattice is given by <fi = 
27r^, where M and N are mutually prime integers. 

For the hexagonal lattice, we can define two triangular sublattices A and B. Let s\, s*2 
and S3 denote the three vectors that connect a site of the sublattice A with its three nearest 
neighbours of the sublattice B, chosen in such a way that s\ A s*2, s*2 A S3 and S3 A s\ are in 
the opposite direction of the magnetic field. We have three hopping amplitudes t^, t\ and 
t\ corresponding to the three directions defined by s*i, s 2 and S3. The index h recalls that 
we have a hexagonal lattice. 

Let us consider a Bloch wavefunction with vector k (0 < k x , k y < -3): 

I ¥ >= 5>**«*4 I >, (2) 

Fi 

where the summation n is over the magnetic unit cell. We only need to consider the above 
region of the momentum, as in the rest of Brillouin zone, the situation can be mapped to 
this case. 

With the Hamiltonian (Hp, the Schrodinger's equation yields the following equation for 
the coefficients corresponding to the energy E h : 

E h u H = t1a1e ie ^ Uri+gi if n e A (3) 

i=l 

E h u H = ]T ^{a^e^-^un-s, if n G B (4) 

i=l 

where a\ = e *. Combining (^) and (0), we obtain the eigenvalue equation for the coeffi- 
cients un of the sublattice A: 
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{{E h f-{{t^ + {t h 2 f + {t } lf)}un = 

iytj = l 

At this point, we define the three vectors Si = s 2 — S3, S2 = s 3 — s\ and S3 = s 1 — s 2 
which connect nearest neighbours of the sublattice A. Using the fact that: 



9ri,n+Si + Qn+Sufi+Si-Sj ~ e ijk ^ + ®n,n+e ijk S k (^) 



with eyfc the Levi-Civita symbol and the flux through elementary hexagons, we have: 

{W-((*f) 2 + (*2 fc ) 2 +(<§) 2 )}«« = 
3 

E WW)"V«** e **^*u^^, (7) 

where a; = e 1 ^. 

Let us then consider the Hofstadter problem on the triangular lattice. With t\,t\,t\ the 
three hopping amplitudes the Schrodinger's equation gives the following equation for the 
coefficients Uft corresponding to the energy E l : 

E'un = £ (fia\e ie *^u H+§i + l&afi^e"*-*-**^) (8) 
i=i 

where a\ = e tk ' Si . 

With the above relations, one can now easily relate the energies of the hexagonal lattice 
to the ones of the triangular lattice: 

(E h f - ((i*) 2 + (i*) 2 + (i*) 2 ) = E* (9) 

where we have the following correspondances of the hopping parameters and of the momenta: 

h \ ~ L 2 L 3 

i.t ihih 

L 2 ~ L 3 L 1 

it ±h±h 

l 3 — t 1 t 2 



(10) 

It is well known that the spectrum of any bipartite lattice is symmetric around zero. This 
can be understood from the fact that the Hamiltonian H can be transformed in —H through 
the following unitary transformation: 

c r if x G A 

(11) 

— c x if x G B. 

The hexagonal lattice is bipartite, while this is not the case for the triangular lattice. There- 
fore the spectrum on the hexagonal lattice is symmetric around zero and this property does 
not hold for the triangular lattice. This result is clearly seen with Eq.(Q). 

The duality relation expressed by Eq.(|9]) and Eq.flTUD allow us to use directly the results of 
Faddeev and Kashaev |J to find the Bethe-Ansatz equations for the hexagonal system. Let 
us recall briefly the results for the triangular lattice. Using a N 3 reducible representation, 
the Hilbert space is a tensor product of three subspaces, each of which has dimension N . A 
wavefunction can thus be written as 

| 4> >=| 4> > ® | 4> >i ® I <P >2 • (12) 

Then, introduce the following notation: 

«*i = e 2 /2 (f 2 C)-V\ai = eTihCY^^l = ef(/ C)-^, 

A = (e 2 f 2 Cf 2 A = {erhC) 1 ' 2 ^ = (e / C) 1/2 . (13) 

One also has 

e* = h-idiCi+i, ft = Ci-xOibi+i, (14) 
with i G Z3. Here C is a constant complex number, C N = (— l)^ -1 ). 
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The eigenvalue equation for the projection Q(p) of an eigenvector | <fi > on the Baxter's 
vector | p >, which corresponds to the point p = (x, £ , £ 1; £ 2 ) on a curve in a four dimensional 
space, reads 

A(x)Q(p) = Q(r_p)A4p) + Q(r + p)A + (p), Q{p) =< <f>\p >, (15) 
where A(x) is related to the eigenvalue E l by 

A(x) = ua aia 2 C + d did 2 + x 2 E f , (16) 
t + and r_ are defined on the coordinates x, £o>£i>£2 by 

T±x = uj ±x ' 2 x, r±& = uT 1/2 &, 1 = 0,1, 2, (17) 

and 

i&Z 3 

= II &( a i d i ~ x 2 biCi)/(£i +1 ai - xbi). (18) 

The zeros pk of Q(p) are given by the following Behte-Ansatz equation 

Q{T-Pk) A + (p fc ) 

Q(r +Pfc ) A_(p fc )" 1 J 

The equations for the hexagonal lattice take the same form except that one has to introduce 
( [TOD into the definitions (pp. 

In the special case where the genus of the curve vanishes and N is odd (N = 2P + 1), one 
may simplify the above equations, so that they can be written out explicitly This special 
case corresponds to the values of momenta: a\ = a\ = a\ = q 1 ^ 2 , with q = uj 1 / 2 . The energy 
for the hexagonal lattice reads: 



(E h ) 2 = (ft) 2 + (t h 2 ) 2 + (^) 2 + (q 1 ' 2 + q- 1/2 )(44 + tffi + t\t 

2P 



2P 



m=l 



+(W(?-?" 1 )(! 1/> -?" 1/S ) E %A (20) 

l<m<n<2P 



where the z\ (1 = 1,..., 2P) are given by the Bethe-Ansatz equations: 

qzi -z m _ _ 1/2 (t\t\z x + q^mtU + q^)(t\t\ Zl + q^) 
m= i 1 i m ^ z, - qz m q (q^H\t\ Zl - l)(qVH h Azi ~ IW' 2 ^ - l\ 1 ] 

In summary, we have developed the Bethe-Ansatz for the Hofstadter problem on the 
hexagonal lattice in this paper. An interesting duality is discovered between the hexagonal 
lattice and its dual partner (triangular lattice). Using this duality relation, we have written 
the Behte-Ansatz equations for the point of the magnetic Brillouin zone where the equations 
of Faddeev and Kashaev take an explicit form. Further work is necessary to obtain more 
information from these equations. 
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by the Swiss National Science Foundation. 
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